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List of Terms; 


F = Forward scattering coefficient for diffuse incidence 
F' = Forward scattering coefficient for parallel incidence 
B = Backward scattering coefficient for diffuse incidence 
B' = Backward scattering coefficient for parallel incidence 
jx = The absorption coefficient for the diffuse flux 
jj.' = The absorption coefficient for the parallel flux 
t = Diffuse flux in the forward direction 
s = Diffuse flux in the backward direction 
I^' = The residue of the primary beam 
J(0) ~ Spectral radiant intensity 
Ajj = Fraction of total absorptance assigned to the diffuse flux 
Ap = Fraction of total absorptance assigned to the parallel flux 
Aj. = Total absorptance 
R = Body reflectance, diffuse incidence 
R' = Body reflectance, parallel incidence 
T = Body transmittance, diffuse incidence 
T' = Body transmittance, parallel incidence 
T ' = Unscattered fraction of the incident beam transmitted 

p 

Tp = Fraction of incident beam diffusely transmitted 
p ' = Unscattered fraction of the incident beam reflected 
Pjj' = Fraction of incident beam diffusely reflected 

= External surface reflection coefficient, parallel incidence 
Tj = Internal surface reflection coefficient 


1 



Introduction 


Beginning with the fact that the units of the absorption coefficient are inverse 
s 

centimeter, it is apparent that some distinction should be made between the ab- 
sorption coefficients /x' and jj. assigned respectively to the parallel and diffuse 
fluxes permeating a medium. It is the purpose of this paper to show that these 
coefficients can be determined simply and directly by goniometric measurements. 
In this approximate treatment, we omit the more general case, in which the ab- 
sorption coefficient may exhibit a directional dependence and consider only the 
mean path associated with the two coefficients. We will assume that the material 
is optically thin, macroscopically homogeneous and that first order scattering 
predominates. Further, we will consider materials in which the scattering is 
concentrated below the critical angle. For transmitting materials these assump- 
tions are not excessively restrictive. It is presumed throughout that the sample 
is optically plane and parallel, that the measured quantity, the radiant intensity 
is monochromatic, and that edge effects are minimal. 

Theory 

The differential equations describing the fluxes permeating a scattering 
medium are given by,^’^-^ 

F' I^' - - Bt + Bs (1) 

B' - Bs + Bt (2) 


dt 


ds 

dx 
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a 

If as predicted, higher scattering orders are negligible, then we can simplify 
(1) and (2) and write. 


dx 


F' I ' 


X 


- /Xt 


(l.a) 


ds 

dx 


B' I ' 


X 


~ /US 


(2.b) 


The attenuation of the primary beam is given by. 


dx 


(m' +B' +F')i; 


( 3 ) 


If we let. 


t 


e-^x F' 





and impose the boundary conditions. 


X = 0 , t = 0 

x = X, t = T'-I' 

X 


the following solution is obtained. 


T' - I ' 

X 


e‘^*F' I 


0 


(a^-A^'-B'-F') 


|g-(Ai'-M+B'+F0x _ 
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which we will approximate by, 


T' - i; = F'X(l-/.x) I,' (4) 

Similarly, from the boundary conditions 

X = 0 , 

X = X , 

we find that, 

R' = B' X(l-/xx)I(,' (5) 

It should be noted that equations (4) and (5) are, as expected, independent of 
n' and that linearity with thickness as required by primary scattering is indi- 
cated for small /^.x. 

Now in traversing a layer dx a fraction of the primary residue is scattered 
and the remainder transmitted, each fraction subject to a particular absorption 
process. In view of this partitioning of fluxes, it is clear from (4) and (5) that 
(F' X + B' X) /L(.x represents the portion of the diffuse flux absorbed. Similarly, 
we will let e"'*'* /x' x , where q' =/x'+B'+F', represent the portion of the un- 
scattered residue absorbed. It is assumed here that is equal to unity. The 
total absor^tevKe is then, 

= (F' X+B' X)/xx + X (6) 

— 

as conservation requires, 

At " 1 - - Pp - Pd (7) 


s = R' 
s = 0 
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The fraction of the total absorptance we can assign to the diffuse flux is obviously, 


(F' X + B'X);ux 


(F' X + B' X) /XX + e”** /x' X 


Substituting from (4) and (5) we have for small /xx, 


(T' - + R') yLiX + /x' X 


(8) 


(9) 


Essentially equations (4) and (5) pertain to a one dimensional medium. 
Consequently, the coefficient, /x is not given in terms of an actual path but rather 
a thickness X.'* It suffices for our purpose if we consider a mean path. The 
probability that a fraction of radiation will be scattered within a solid angle 


du> - s in (9 d(9 dcp 

in the direction 6 is given by,® 


S(0, «p) da> 

where S(9, cp) is the scattering function. Therefore, for an axial symmetric 
case the mean direction is, 


/V f 0 S(0) doj 

<e) = J V ~ (i<» 

Since the scattering function is directly proportional to the radiant intensity, 1(9) 
for either conservative or non-conservative scattering it suffices for most cases 


5 



to plot the flux. 


J (^i) sin 6'. A6* 

as a function of 6 and determine the mean by inspection. For example in figure 
(3) the indicated mean of about 28.5° is in good agreement with the calculated 
value. It follows then that if (i9) is the mean direction, then X/^cos 6 ^ is the 
mean path. Hence, 



jjX - X 

(11) 

Consequently, (9) can be written as. 


A — 

[(T'-i;)+R']/<cose> 



[(T'-I,')+R']/<co,e> ti; ^ 

(9.a) 

If we now consider the boundary effects resulting from the refractive dis- 

continuity, the parallel and diffuse fractions of radiation transmitted and 

re- 

fleeted are, 



T ' = 

P 

(l-ro)2l/(l-r,2i'2) 

(12) 

^d' " 

(l- ro)(w +VFo I'2) 

(13) 

^p' = 

To + 

(14) 

Pd 

(l- r„) (v +wto I')/(l- r<2 I' 2 j 

(15) 
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where, 


(r-I')(l-r,RUr,R'T 

w - — ; 

^ (l-r^ R)2- r/ T2 

R'(l-r2R)+r2(T'-I')T 

For transmitting materials, r^, r^, R and R' are often relatively small. We can 
therefore rewrite the previous equations as. 


^p' (l- ‘‘o)'* I' 

^d' ^ (i~ ^o) (l - 

Pp ro+ >'o (l 

Pd (l- ^o) (l - 

Now assuming that. 



(12.a) 

2)(T'-I') 

(13.a) 

)2 I'2 

(14.a) 

2)»' 

(15.a) 


~ /XX + /Lt' X 

and substituting (12.a), (13.a), and (15.a) in (9.a) we have approximately, 


(■^D +Pd')/^^os 0 ) 
i'^D + Pd )/ (cos 0) + Tp' ^ 

Similarly we can show that, 

^ (V + ^d)/(cos 0> +Tp' ^ 


(16) 


(17) 
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Procedure: 


Since the measured flux has been refracted and reflected at the boundary 
and our immediate concern is with the internal angular distribution, it is neces- 
sary that we determine in terms of the measured radiant intensity If 

for the moment we ignore reflection losses and consider only the refractive 
effect at the boundary, then the flux Incident must equal the flux emerging at the 
boundary. Therefore, 


or 




J ( 61 ) _ doj 

J(6) dc<j' 


( 18 ) 


The ratio of intensities is then equal to the inverse ratio of their respective solid 
angles. From Snell's law 


n sin 0 - n' sin 0 ' 


we find that. 


n cos 0 d0 - n ' cos 0' d0' 


Hence, 


doj' _ sin 0 ' d^ ' dcp _ j cos 0 
do; ” sin0d6'dcp “ cos0' 


which can also be expressed as. 


do) 


n^ /l - sin^ 0'/ n^ 
cos 0 ' 


( 19 ) 
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For normal incidence this reduces to the more familiar relationship,® 


da' - da 

The two intensities are related by, 

fn^ ( 1 - sin^ 0 Vn^ ) , 

^ cos 

If we now include reflection losses at the boundary, 


(S) 


r r> a (l - sin^ '/n^ ) 


( 20 ) 


where B.(6) is the general Fresnel expression. 


R(d) 


1 sin^ ( 6 ' - (9) J cos^ (6 


2 • 2 


1 + 


sin (0 + 6' ) cos^ (6 




( 21 ) 


The procedure is now rather straightforward. From the measured intensity 
J(g), J( 5 )is calculated using equation (20). As an example, the two intensities 
are shown in Figure (1) for a 3.8 mm thick Irtram I sample at 1/i.. The Fresnel 
reflection coefficient for a refractive index of 1.38 is plotted in Figure (2). Now 
as shown in Figure (3), the flux da is plotted and a mean direction determined. 
Returning to equation (7) the total absorptance can be deduced by integrating the 
forward and back scattered flux, ^ 


= 277 




(^0 


sin 


d. A0 


i=i 


I o' 


( 22 ) 
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and 


7T Ln . sin (9 . IS6 

^ L r? 


i=7T/2 


(23) 


Alternatively, t^' and can be determined directly with an ellipsoidal or inte- 
grating sphere photometer. The transmittance, r^' is readily determined. The 
’’eflectance, however, is not generally susceptible to direct measurement. 

We can, however, solve equation (12. a) for I' and then substitute in equation (14. a). 
So that, 

Pp " + ((24) 

The rest is a matter of substituting in the appropriate equation. 

Results : 

We can in some measure validate both the assumptions and method outlined 
here for determining the absorption coefficients p and p' . We begin by restating 
equations (13. a) and (15. a), 

(1-n) " (Ti;y (25) 

R' = PD'/(l-'o){l-n) ( 26 ) 

Again, following Ryde we note that as Q - 1 andP - 0, the equations for the body 
transmittance. 


QK + sinhKX ^ 

T - (yLt +B) sinhKX + K cosh KX ■ (Q 


( 27 ) 
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and reflectance, 


Pe"'’'’‘K + QBsinhKX 
(m + B) sinh KX + K cosh KX " ^ 


(28) 


can be written as, 


and 


T' - T 


K 

(m +B) sinhKX + KcoshKX 


(27.a) 


R' - R 


B sinh KX 

(fi. +B) sinhKX + K cosh KX 


(28.a) 


where, 


KX = /mX (ijX + 2BX) 


(29) 


Significantly, both R and T are a function of and B. Furthermore, an approx- 
imate value of the backscattering coefficient, B is obtained from the ratio, 


Expanding sinh KX, 


R BX sinhKX 
T “ KX 


R 

T 


B 

K 



(KX)3 (KX)5 



we find that. 


R 

T 


BX 
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I 


oi' substituting from (25) and (26), 


BX 








( 30 ) 


The following measured values were obtained, r^' = 0.23, = 0.45, 

Pd' = 0.0725 and from (24) p^' = 0.0274. Assuming = 0.15 and given that 
Fp = 0.026, (cos 6’} = 0.879, we find from (30), (29) and (16) that BX = 0.1114, 

KX = 0.246 and px = 0.1586. Inserting the above values in (25) and (27.a) we 
find that T = 0.785 and T = 0.769, indicating agreement within 2%. Similarly, 
from (26) and (28. a) we find that R = 0.08624) and R = 0.087571, indicating agree- 
ment within 1.5%. The measured values for other thickness of Irtram I further 

Q 

substantiates the values obtained for BX and px reported here. 

Conclusions : 

An approximate method has been outlined for determining the absorption 
coefficient p and p ' which appears to be applicable to optically thin materials 
where primary scattering predominates. The method has at least the advantage 
of following directly from a fairly conventional measurement. The coefficients 
are given in terms of the reflected and transmitted fluxes weighted by their 
mean paths. Some of the assumptions made, particularly a method for deter- 
mining the interval reflection coefficient, will be the subject of another paper 
now in preparation. 
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